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Abstract 

The production and decay rate of massive sterile neutrinos at finite temperature 
receives next-to-leading order corrections from the gauge interactions of lepton and 
Higgs doublets. Using the Closed-Time-Path approach, we demonstrate that the 
perturbatively obtained inclusive rate is finite. For this purpose, we show that soft, 
collinear and Bose divergences cancel when adding the tree-level rates from 1 -f-)- 3 
and 2 -f-)- 2 processes to vertex and wave- function corrections to 1 -H- 2 processes. 
These results hold for a general momentum of the sterile neutrino with respect to 
the plasma frame. Moreover, they do not rely on non-relativistic approximations, 
such that the full quantum-statistical effects are accounted for to the given order 
in perturbation theory. While the neutrino production rate is of relevance for 
Leptogenesis, the proposed methods may as well be suitable for application to a 
more general class of relativistic transport phenomena. 
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1 Introduction 



Scattering processes between particles at finite temperature play an important role in 
Early Universe Cosmology in view of applications such as Baryogenesis and the produc- 
tion of Dark Matter. Besides, certain properties of hot and dense strongly interacting 
matter, such as the photon production rate and transport coefficients rely on the knowl- 
edge of these scattering rates. 

It turns out useful to distinguish between the scatterings of four massless particles, 
which are at leading order restricted to 2 f-)- 2 processes and the case where one of these 
particles is massive, which opens up the phase space for 1 <f->- 2 and 1 O 3 processes as 
well. In the case of the scattering of massless particles only, an important complication 
can arise from the t-channel exchange of a massless particle, leading to a Coulomb 
divergence at tree level. This divergence is mitigated by the screening in the plasma, 
and technically, it can be removed by a resummation of thermal self-energies within 
the propagator of the mediating particle [l|-[3[. When one of the scattering partners 
is massive, this divergence no longer occurs, but instead, there are soft and collinear 
divergences. 

The soft and collinear divergences are of the same type as e.g. those familiar from 
QCD corrections to the quark pair production from an off-shell photon. While the 
production rate for quark pairs and a soft or collinear gluon becomes non-perturbative 
below a certain transversal momentum scale between the gluon and the emitting quark, 
the inclusive cross section for pair production with and without gluons is perturbatively 
well defined, due to the cancellation of soft and collinear divergences between tree-level 
and loop diagrams. In the following, we refer to the soft and collinear divergences 
collectively as infrared (IR) divergences. 

It is natural to consider the same problem at finite temperature. Recently, this 
matter has received attention in the context of Leptogenesis The quantity of 

interest is the relaxation rate of sterile right-handed neutrinos toward equilibrium 
from decays and inverse decays of Higgs bosons and leptons. When the mass M of is 
small compared to the temperature T, radiative corrections are of leading importance, 
because the 1^2 processes are kinematically suppressed. In the opposite case, M 3> T, 
radiative corrections are subdominant to the leading 1^2 rates, but yet, it may be 
interesting to accurately know the size of the next-to-leading order (NLO) corrections. 
Furthermore, the situation M S> T corresponds to strong washout, which is in type-I 
see-saw scenarios perhaps favoured in the light of the observed mass-scale of the active 
neutrinos. Whether the 1 -H- 2 processes are kinematically suppressed due to M <C T 
or not, we refer to the radiative corrections from Standard Model gauge interactions as 
NLO in this paper. The problem of the NLO corrections to production for M ^ T 
has recently been resolved js], lij , making use of the fact that Maxwell statistics is a good 
approximation in the non-relativistic regime and that one may assume that the scale of 
the three momentum of is much below the temperature T. Related to this matter 
are also works on the Drell-Yan production of dileptons from a strongly interac ting 
plasma [7-10|, and some useful calculational techniques are developed in Ref. 11, 12 
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where the spectral function of hot Yang-Mills theory is studied. 

However, it still remains interesting to demonstrate the cancellation of IR divergences 
between tree- level and loop corrections at NLO, keeping full account of relativistic effects 
and of quantum statistics. This is the main result of this work. The proof of the 
cancellation of the IR divergences that is presented here relies partly on analytically 
isolating them from the integrals (for the wave-function corrections) and partly on a 
rearrangement of the integrand such that only integrable singularities remain (for the 
vertex corrections). Therefore, this proof is already suggestive of a method of numerically 
calculating the production rate, a task that we will pursue in a separate work. 

In order to compute the relaxation rate for N at finite temperature, we choose to 
use the Closed-Time-Path (CTP) method 13|, In this formulation, Schwinger-Dyson 



;hXCT 

equations can be derived |l5l-ll7j| (a particular subset of which are known as the Kadanoff- 
Baym equations) that describe the time evolution of the Quantum-Field-Theoretical 
non-equilibrium system. For this reason, the CTP approach has recently been applied 
to Leptogenesis in different parametric regimes, and various new effects and corrections 
have been derived 18 -S^] • In the CTP formulation, the relaxation rate is directly 



related to the spectral (or anti-hermitian) self energy. As we assume that the fields 
that participate in gauge interactions are in thermal equilibrium (which applies to all 
fields within the loop diagrams that contribute to the self energy of N in the present 
approximations), we notice that the relaxation rate can be obtained from equilibrium 
field theory as well. Here, we choose the CTP approach because of its connection to 
non-equilibrium field theory formulations for Leptogenesis. Moreover, the relation of the 
spectral self energy to the relaxation rate can be easily inferred from the kinetic equations 
that are derived from the Schwinger-Dyson equations, (c/. Ref. 2^), and finally, the 



CTP approach provides very efficient Feynman rules that are perhaps easier to apply 
than certain finite-temperature cutting rules used in Ref. [l]. 

The methods that we employ in order to demonstrate the cancellation of IR di- 
vergences are very explicit, which has the disadvantage that they rely on a number of 
technical details on the evaluation of phase-space and loop integrals. On the other hand, 
the explicit demonstration of the cancellation readily suggests a method of practical 
evaluation of the neutrino relaxation rate to be pursued in a future work. In order to 
communicate the main points of this work to a reader who is not interested in the tech- 
nical details and to give a first overview to a reader who is at least partly interested in 
these details and their reproduction, we provide in Section [2] a summary of our method 
and of the cancellation of the IR divergences, that gives reference to the main results 
that are worked out in the technical Sections. We begin Section [2] with a part containing 
general prerequisites and remarks about the CTP approach and the neutrino relaxation 
rate ( Section [2TT!) . As suggested by the different diagrammatic contributions to the relax- 
ation rate, which also lead to different calculational methods, we further present separate 
parts about wave-function type (Section 12.21) and vertex-type (Section 12. 3p corrections. 
In Section [31 we present the main results for the IR- and ultraviolet-finite neutrino pro- 
duction rate, and the the full calculational details are then presented in Section H] for the 
wave-function corrections and in Section [5] for the vertex corrections. We conclude with 
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Section [61 



2 Summary of the Calculation and Cancellation of 
Collinear, Soft and Bose Divergences 

2.1 Kinetic Equations, Spectral Self Energy and the Relaxation 
rate for Right-Handed Neutrinos 

We are interested in the production rate of right-handed neutrinos A^, or more precisely, 
in their relaxation rate toward thermal equilibrium. The self energy for N is given by 
t^]\f. At one- loop order, it originates from Yukawa interactions Y (that we choose here 
to be real) with lepton doublets i and Higgs-boson doublets (j). Besides, receives 
higher order corrections due to the SU(2) and U(l) gauge interactions in the symmetric 
Electroweak phase at high temperatures. In particular, the model we consider is given 
by the Lagrangian 

C = CsM + ^^jv(0 - M)^jv - r^,0tpj^^^ _ Yt/j^^P^ij, , (1) 

where £sm is the Lagrangian of the Standard Model, ipN/ are the spinors associated with 
and i and Pl,r = (1 T7^)/2. (Group indices of SU(2) and antisymmetric tensors that 
may be necessary to form invariant products are suppressed throughout this paper.) On 
the spinor associated with A^, we impose the Majorana condition i/)^ = Cipjf = ipN, 
where C is the charge-conjugation matrix. The gauge couplings are given by by g2 
for SU(2) and gi for U(l), and moreover, we denote by (7^ = 2 the dimension of the 
fundamental doublet representation of SU(2). As the Higgs doublet and the lepton 
doublet have the same Electroweak charges, it is useful to define the recurring factor 

G = l9l + \9!. (2) 

Besides, there are sizeable corrections from top-quark loops and self-interactions of the 
Higgs boson. While we do not include the latter explicitly, their treatment should be 
straightforward, given the discussion of the wave-function corrections in Sections 12.21 
andU 

Since we assume the leptons, Higgs bosons and gauge bosons to be in thermal equi- 
librium, the relaxation term in the kinetic equation for the right handed neutrinos is 



proportional to their deviation from equilibrium |24i |: 

|/iv(k) =\ I ^sign(fco)tr [i^>{k)iS<{k) - i?< (A;)i^>(fc)] (3) 



1 f dh 







^ , ^ ^sign(fco)tr mi{k) 27r5(fc2 - m^)5/Ar(k) 
Z I Ztx L J 



2F'^ 



Nik) 5/^(k), 
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Figure 1: The self energy [^jy]^*^ 



where = ±a/I?+M^ in the last term and where SJnO^) is the deviation of the 
distribution function /Ar(k) of right-handed neutrinos from the Fermi-Dirac distribution. 
We assume spatially homogeneous conditions and work in Wigner space, where the two- 
point functions are Fourier transforms with respect to the relative coordinate. Therefore, 
there remains a time dependence, but we however suppress the explicit time coordinate 
in the arguments of the Wigner functions and distribution functions. As we calculate the 
self energies in this paper assuming thermal equilibrium for the internal propagators, we 
make the remark that the approach we use here is effectively identical to the real-time 
formulation of equilibrium field-theory. In Eq. (I3D, we have made use of the definition 
of the spectral self-energy (lAldl) and of the KMS relation for equilibrium two-point 
functions (lA4p . Notice that when ^"^{k) is known in equilibrium, the KMS relations 
immediately yield t^'^ik) and therefore the relaxation rate. 

In order to calculate ^jv(^) when the mass M of is non-zero, we perform a loop 
expansion. For this expansion to converge, it is crucial that we calculate an inclusive 
production rate with or without the production or absorption of a real gauge boson, such 
that IR-divergent contributions from real emissions (2^2 scatterings and 1^3 decays 
and inverse decays) and loop corrections (to the 1 ■v^ 2 decays and inverse decays) cancel. 
Note however that the perturbation expansion breaks down in the limit M — )■ 0, because 
of the t-channel exchange of massless fermions [11-0]. We will deal with this matter in a 



separate publication [35 



The Yukawa coupling between the right-handed neutrino A^, the lepton doublet i and 
the Higgs doublet gives rise to the leading-order contribution to the self energy of the 
right-handed neutrino, that is represented diagrammatically in Figure [1], 

ittip)]'^'' = 9^Y' J ^(i5f(p + A:)iA;''(-A:) + C[i5f(-p-A;)]*CtiAj«(A;)) , (4) 



where a, 6 = ± are CTP indices (we follow the conventions used in Ref. [16[) and t stands 
for a transposition of Dirac matrices. (We indicate the various contributions to t^j^ by 
superscripts on square brackets, for the sake of readability. Notice that we deviate from 



the conventions of Ref. [16|] in the detail that we define the fermionic self-energies as 
^, while in Ref. 16] the same quantities are defined without a slashj^]) The tree- level 

^Here it is understood that ?5 may involve chiral projection matrices -Pl.r, for instance for left-chiral 
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propagators for the lepton and the Higgs fields are given by Eqs. (lASp and (IA6p . Exphcit 

V . A 1 LO 

expressions for 



24 



can be found in Ref. 

The main purpose of this work is to provide a method for calculating the NLO 
contributions. These may be categorised in two-particle-reducible wave-function type 
corrections (superscript WV, Figure [3]) and two-particle-irreducible vertex-type correc- 
tions (superscript VERT, Figure H]). Due to the t-channel divergences, it turns out that 
the two-particle- reducible contributions diverge for M — )■ 0, and instead, a two- particle- 
irreducible one- loop diagram with resummed propagators should be calculated |35|. 

In all our calculations, we denote by 



= (M, 0,0,0) 



(5) 



the momentum of the right-handed neutrino in the Centre of Mass System (CMS). The 
relative motion with respect to the plasma is accounted for by a plasma vector 



(6) 



where p° and p are energy and three-momentum of the right-handed neutrino in the 
plasma frame (where the plasma is at rest). Moreover we note that 



(7) 

which follows from the relation (lASp . such that our choice for p with p° > does not 
imply a restriction of the results. 



2.2 Summary of Cancellation of IR Divergences in the Wave- 
Function Diagrams 

The wave-function type contributions can be written as 



1 WV 



iSi'>{p + k)iA>{-k) + \S>{p + k)\A'r{-k) 



(i)>/ 



(8) 



the self-energies 



iwv 



where iSf^ and iA^^'' are the one-loop corrections to the propagators, i.e. 
with external legs (c/. Figure [2]). The particular diagrams contributing to [^jyj 
are therefore two-particle reducible, as apparent from the diagrammatic expression for 
in Figure [31 It would be a simple matter to use fully resummed propagators 
rather than the single-loop insertions, however it is interesting to verify that perturbation 
theory is well defined without resorting to a resummation. 



fcrmions ^ = Pr7^S^. 
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Figure 2: The amputated diagrams in this Figure represent the gauge contributions to 
't^i (first diagram) and (the sum of the second and third diagram). When these 
diagrams are not amputated, they correspond to iS^l'^ (first) and iA^^^ (sum of second 
and third). We refer to the first and second diagram as sunset diagrams, to the third 
as seagull diagram. Due to the CTP Feynman-rules, the seagull diagram contributes to 
nf , but not to nf . 



The rate can be decomposed into a contribution T from gauge-boson radiation from 
the lepton and a contribution B from gauge-boson radiation from the Higgs boson, such 
that 



tr 



Nip) 



1 wv 



(M 



tr 



rtr 



J^ + B, 

iS^'>{p + k)iA>{-k) - i5«<(p + k)iA<{-k) 



iS>{p + k)iAll>{-k) - iS<{p + k)iAl"-{-k) 



(i)</ 



(9a) 

(9b) 

(9c) 



The diagrams that represent [l^^v]^^ ^ind its decomposition in J-" and B are shown in 
Figure O 

In order to explain how the cancellation of IR divergences works, we focus on emis- 
sions of gauge radiation from the Higgs boson, that is captured by B. Radiation from the 
lepton can be treated along the same lines, but is technically slightly more involved due 
to the spinor structure. In Section [5l we present all the necessary details on gauge radi- 
ation from the lepton. The one-loop correction to the scalar propagator can be written 
as 



iA«<'>(fc) =2Af'\k)f<'>{k.u), 
f<{k ■ u) =Uik ■ u) , f>{k ■u) = l + U{k ■ u) 



(10a) 
(10b) 



where A'^^^\k) is the one-loop correction to the spectral function, cf. Figure [T] and 
Eqs. (lAip . As we assume that thermal equilibrium is forced by the gauge interactions, 
throughout this paper, we take for the distributions of Higgs bosons and gauge bosons 
ftf) and Ja the Bose-Einstein distribution and for the distribution of leptons fe the Fermi- 
Dirac distribution. Nonetheless, we keep the subscripts A, cf) and £ in order to indicate 
the origin of the individual distributions. 
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+ 





Figure 3: The self-energy term [^^]^^. The first and second terms constitute the 
contribution J^, the third to sixth the contribution B to tr[|*^'^]^^. 



Now in thermal equilibrium, the one-loop correction to the spectral function can be 
expressed as |36|, [37 



1 



a;^^^(A:) =- {[A^{k)m^ + uf] - [A^(fc)]^[inf - n^]) 

=1 {[A^ik)r - [A^{k)r) in? + i {[A^{k)r + [Ai{k)r) 



(11) 



where 11^ ' are the hermitian and spectral self-energies of the Higgs boson with a gauge 
boson in the loop. [Cf. Eqs. (lAip for the definition of spectral and hermitian two-point 
functions.] In principle, it is important to assume here thermal equilibrium, because oth- 
erwise, there would also occur the product of a retarded and an advanced propagator, i. e. 
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a pinch singularity, which leads to an ill-defined integral because of a double pole above 
and below the real axis. However, it is demonstrated in Ref. j36| that a resummation of 
all loop insertions (in contrast to the single loop insertion performed here) removes the 
pinch singularities and leaves behind only the equilibrium contributions to the propaga- 
tor. In Ref. |37|], it is shown that when performing the gradient expansion, that occurs in 
the Wigner space formulation, to all orders, the non-equilibrium part of the propagator 
takes the intuitively expected resummed form as well. Recently, in Ref. 38|], it has been 
proposed instead to deal with this issue by using a two-momentum representation of the 
two-point function in contrast to the single-momentum representation in Wigner space. 
The retarded and advanced propagators that appear in Eq. f ill I) are given by 



R,A 



iA 



such that in the distributional sense, 



k"^ ± i sign/c'^e 



i^ik)? - [A^{k)r = - 27ri 5^ 5(A:° ± |k|) 



1 d 



(12) 



(13a) 



(13b) 



Using above identities together with the expressions for the tree- level propagators (]A6fA7p . 
the one-loop corrections due to real and virtual gauge bosons attaching to the Higgs bo- 
son are 



d^k 

(2^0^ 

d'^k 
(2^ 

d'^k 



tr yfis>{p + k)i/^f''{-k)-fis<{p + k)i/s.'^;''^{-k) 

J2 /2AJ)-^(A;) [feiip + k) ■ u) + U{k ■ u)] 

fcO = _pO-|-|k| 



(1)</ 



(14) 



(27r 



E 



P 



27r^5(A;° ±' |k|) 



fcO = _pO-|-|k| 

X [h{{v + k)-u) + U{k-u)] 



1 



1 d 



Ak^"^ dk^ Ak"^'^ 



K(k) 



In order to identify the IR-divergent contributions, it is useful to split the self-energy 
into parts that are present for T = with a superscript vac and parts that vanish for 
T = with a superscript T 7^ 0: 



(15) 



The self-energies II^^ and II^^ give rise to two different contributions to B. To those 
from n^'^, we refer to as wave- function contributions, and to those from 11^ as scattering 
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contributions. Notice that the scattering contributions include besides the tree-level 
2 o 2 scattering rates the 1 o 3 decays and inverse decays as well. 

We now briefly discuss the various contributions to B resulting from this splitting 
and explain how the IR divergences present in the particular contributions can be seen 
to cancel eventually: 

• Finite temperature contributions: The contributions from 11^ ' to Eq. flMl) lead 
to IR-finite integrals. It is well known that the hermitian part of the thermal self- 
energy 11^'^^^ is finite, cf. the explicit expression f l27l) below. Therefore, it gives 

rise to a finite contribution to the integral (JH]). For the contribution from 11^'"^^°, 

we show that we can express Uf^^\k) = k^h{k), where h{k) is a continuous 
function, cf. Eqs. ( 136|) below. Therefore, the first order singularity at |k| = k^ (or, 
equivalently, |k| = is integrable in the principal value sense. 

• Vacuum wave- function contributions : It is well-known that the vacuum self-energy 
Yj^H,vac ^ scalar field with a gauge boson and a scalar boson in the loop is IR- 
divergent. This divergence can be regulated by a fictitious gauge-boson mass A. 
More precisely, the A dependence is oc log A and given by Eqs. fl24b|25l) below. 

• Real gauge boson emission (scatterings) in the vacuum : On the other hand, 11^'™'^ (A;) oc 
fc^, such that the resulting singularity at |k| = k^ is again first order. Now however, 
this gives an IR-divergent contribution to the integral, because n^'™'^(fc)/fc^ is not 
continuous for k'^ = 0. Therefore, we introduce for this term the gauge boson mass 

A as well. The contribution from n^'™^ to B adds to the rate of 1 — )■ 3 particle 
decays of A^. Because n^(fc)/fc° itself corresponds to a decay rate of a scalar of 
mass-square fc^ into a massless scalar and a gauge boson of mass A, there is a 
kinematic threshold that can be expressed in terms of a Heaviside "(^-function, cf. 
Eq. ( l28l) below. As a result, this contribution to the integral has an IR divergence 
oc log A as well. When we perform an integration by parts, we can isolate the IR 
divergence and compare it with the one from wave-function corrections. 

• Cancellation of IR divergences from wave-function and scattering corrections: Once 
the IR divergences are isolated in terms oc log A by performing the integral over 
the 6 function in Eq. (fT4|) (for the virtual contributions from 11^'™'^) or through 

integration by parts (for the real contributions from 11^'™'^), one can see that the 
dependences on log A cancel in the total result. For this purpose, compare Eqs. fl33|) 
and (|38|) . This can be viewed as a consequence of the fact that 

- A')^^i^nf--(fc) = sign{k')Uf-{k) , (16) 

which is no accident, because n^'™'^ and n^'™*^ are real and imaginary part of the 
same analytic self-energy. 
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The contributions from gauge-boson radiation from the lepton are discussed in detail 
in Section 14.21 For the vacuum contributions to the spectral and hermitian self-energies, 
the discussion goes along the same lines as for the radiation from the Higgs boson. In 
addition, there are IR divergences when the thermal part of the hermitian self-energy is 
evaluated on shell, which are matched by IR-divergent contributions from the thermal 
part of the spectral self-energy. We note here that these thermal contributions which 
lead to IR divergences can be identified with the self energies that one obtains in the 
hard thermal loop (HTL) approximation. Otherwise, the method of demonstrating the 
cancellation of these divergences is very similar to the one applied to the cancellation of 
the IR-divergent vacuum-contributions. 

2.3 Summary of Cancellation of IR Divergences in the Vertex 
Diagram 

The method relying on the relation flTBl) between the spectral and hermitian self-energies 
that we employ for the two-particle-reducible wave-function contributions obviously can- 
not be applied to the present case of the vertex diagram. Figure HI One may again use 
a gauge-boson mass in order to regulate the soft and collinear divergences. Due to the 
multiple integration boundaries that appear for virtual corrections and real emissions, 
one yet faces the problem of matching the particular contributions that lead to a cancel- 
lation of the IR divergences, i.e. of the dependence on the fictitious gauge-boson mass. 
This task can be facilitated by transforming the integrals and by arranging and adding 
the particular integrands in such a manner that there remains a single integrand that 
manifestly evaluates to a finite answer. This is the method pursued in the present work. 
One could then either go back to evaluate the particular divergent contributions with a 
regulating gauge-boson mass, since it is then clear that the dependence on the regulator 
eventually cancels from the finite result. Alternatively, one may directly integrate the 
total integrand without introducing an IR regulator, what again leads to a finite result. 

In order to accomplish the task of arranging the particular contributions to the vertex- 
type self-energy into a manifestly IR-finite integral, we focus on the reduced Wightman 



handed neutrinos tr[^i|5^(p)]/(4p°) (for j9° = a/p^ + M^) [cf. Eq ([3]) and the defini- 
tions ( lAll) ]. The thermal decay rate tr[^i^^(p)]/(4p°) may be directly inferred from the 
KMS relation (jAil), i^^(p) = -e'^P-"i^^(p), and the relaxation rate tr[^|i^(p)]/(2p°) 
can then be found using the definitions (lAlll . Besides, we assume first that p° = M > 0, 
while the solution in the case p° < can be obtained from relation (JTj). 

The CTP expression for the 2-particle-irreducible vertex-type correction to the right- 
handed neutrino self-energy, represented graphically in Figure HI is 




This quantity is proportional to the production rate of right- 



'N 



n VERT 



(p) =2g^GY^ 



I 



jJ2cdiSr{k + q)riAfM (17) 



(2vr)4y (27r) 



c,d 



x[2ik-p) + qriSfik)l 



{-k + p)iA^{-k + p-q) , 
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+ 




Figure 4: The self energy term 



IVERT 





Figure 5: Vertex ([^jv]™'*! 



and scattering ([|!ljv]'^'^'^) contributions to [t'N] 



IVERT 



where the tree-level propagators are given by Eqs. (1A5IIA6|A7I) . 

The reduced correction to the Wightman self-energy tr[|*i^^(p)]^™'^ follows from 
the a,b = +, — contribution. Furthermore it is useful to shift in the two terms where 
c = — the momenta as — )■ + g and q — )■ —q, to make use of the cyclicity of the trace 
and of the definitions flAlcfAlep . such that one obtains 



VERT 



d^k d^q 



[2{k - p) + q^ti 



(18) 



(27r)4 (27r)4' 

X \\Sf{k)Y-iS<{k + g)iA<(-A: + p)iA>,(g)iAj {-k + p- q) 

+iSf{k)YiS<{k + g)iA^(-fc + p- g)iA J {-k + p)iA^,M 
+iSf{k)YiS<{k + q)iA<{-k + p- q)iA^{-k + p)iA^M 
+iS[{k)riS< {k + q)\A<{-k +p- q)iA^ i-k + p)iA^M 

Notice that the Wightman and statistic propagators (superscripts <,>,F) are purely 
on shell, as they are proportional to an on-shell 6 function, whereas the hermitian prop- 
agators (superscript H) correspond to principal values. It is therefore easy to see that 
in this decomposition, the four terms in the integrand correspond to the following rates: 

• The first term describes the interference of two tree-level 2^2 scattering or 
1^3 decay diagrams, where in one of these, the gauge boson radiates from the 
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lepton, in the other one from the Higgs boson. We denote these contributions 
collectively as scatterings (even though for 1^3 processes this term does not 
strictly apply) with the superscript sea, tr[^i^^(p)]'^'^'^. Diagrammatically, this 
contribution corresponds to a cut through the Wightman propagators [those with 
superscripts <, >, cf. Eq. flA2p ]. which are purely on shell, cf. Figure O 

• The second term yields a correction to the 1 4-)- 2 process. It is the interference 
of the part of the vertex correction where the gauge boson is on shell with the 
tree-level amplitude. We denote this contribution by tr[^i^^(p)]™''*. 

• The third term yields a correction to the 1 2 process as well. It is the interference 
of the part of the vertex correction where the Higgs boson is on shell with the tree- 
level amplitude. We denote this contribution by tr[|ii|!l^(p)]2°'^*. 

• The fourth term yields a correction to the 10 2 process as well. It is the inter- 
ference of the part of the vertex correction where the lepton is on shell with the 
tree-level amplitude. We denote this contribution by tr[^i|^^(p)]™''*. 

The vertex corrections can be represented by cuts through the Wightman propagators as 
well, cf. Figure [51 We furthermore emphasise that all of these individual contributions 
contain the full quantum statistical factors for the on-shell particles. 

For all of these four individual contributions, IR divergences can occur whenever the 
hermitian propagators go on shell. In order to identify the location of the IR-divergent 
contributions to tr[|ii^jy(p)]*^'^'^, it is therefore useful to parametrise the integrals by vari- 
ables that are proportional to the invariant momentum squares of these propagators. 
Here, we choose x = {—k + p — q)'^/M and y = k'^/M with the momenta as parametrised 
in Eq. f|T8l) . [Cf. Eqs. (177]) for the parametrisation in terms of the momenta chosen in Fig- 
ure O] In addition, as we are working in a homogeneous finite-temperature background, 
there eventually remain two non-trivial angular integrations. The kinematic constraints 
that are forced within the CTP approach by the on-shell delta functions within the 
Wightman and statistic propagators then imply that the integrand for tT[^it^^{p)Y^^ has 
a limited support, which is illustrated in Figure [61 Notice that the resulting sub-domains 
have simple kinematic interpretations that are given in the Figure caption. The lines 
where x = or ?/ = correspond to the location of collinear divergences, while the point 
where x = y = to the soft divergence, where the emission of a zero-momentum gauge 
boson leads to a Sudakov logarithm square. Because the collinear fringes where x = 
and y = are also the boundaries of the integration domain for the real emissions, the 
singularities from the hermitian propagators cannot be integrated in the principal value 
sense but lead to logarithmic IR divergences. In addition, there are Bose singularities 
from the gauge boson distribution at x = y = and from the Higgs boson distribution 
at(x,i/) = (0,M/2). 

These divergences cancel with IR divergences from the corrections tT[pit^j^{p)]i'^2%- 
What we demonstrate in this work is that the integrals defining tr||ii^^(p)]^'^2*3 can be 
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Ill 




III 



Figure 6: Regions of integration for tT[pit^j^{p)Y'^^, corresponding to tree-level 1^3 
and 2^2 processes. The coUinear singularities coincide with the solid blue boundaries, 
X = or y = 0. Region I corresponds to 1 -H- 3 processes, Region II to the interference 
of two t channel amplitudes, the lower branch of Region III to the interference of one 
amplitude with a lepton in the s channel with one with a Higgs boson in the t channel 
and the upper branch of Region III to the interference of one amplitude with a lepton 
in the t channel and one with a Higgs boson in the s channel. 

transformed and added to tr[|ii^^(p)]^"^^ in such a manner, that the resulting total inte- 
grand in tr[|*i|i^(p)]^^^'^ has only integrable singularities after all. This procedure is also 
suggestive of a practical method for numerical integrations or analytical approximations, 
that can be described as follows: 
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• Parametrise the integral tT\pit^j^{p)f'^^ in terms of x,y as described above and as 
illustrated in Figure |6] (Section l5.ip . 

• Parametrise the integrals tr[|Si|^^(p)]™|*3 by x and y as well, such that these corre- 
spond to 1 /M times the momentum squares of the Higgs and the lepton propaga- 
tors in the sub-diagrams that correspond to virtual vertex corrections. Add various 
contributions such that the integration area for the virtual corrections covers the 
quadrant x > 0, y > (Sections 15.21 15.31 and 15.41) . 

• For the real corrections, fold Regions II and III (either or both x,y < 0) to the 
quadrant where x,y > 0. 

• Now add all contributions. The resulting integrand J'^°^^^{x,y), Eq. (I116bp . only 
contains integrable singularities (Section l5.5p . 

• The integral is not yet convergent for large values of x, y, what corresponds to 
an ultraviolet (UV) divergence. We obtain a convergent integral by subtracting 
a term that corresponds to the vertex correction to the vacuum decays N ^ 
weighted by the thermal distributions for i and 0, jf^'^^'t'^'^^ Eq. f ll30ap . As this 
contribution is IR divergent, we cancel it by adding a correspondingly weighted 
rate for 1 — )■ 3 decays, jf'^ca.vac^ (11321) . The subtracted contributions must be 
added again to the final results, but for these particular terms, we can perform 
the integration over dxdy analytically and thus isolate and renormalise the UV 
divergence (Section l5.6p . 

The point stating the presence of integrable singularities only requires a proof, that 
we present in Section [5] and that at least in parts elucidates how the cancellation of IR 
divergences works in the finite-temperature background. The various contributions to the 
integrand obviously factorise into kinematic and quantum statistical terms. We perform 
an expansion of the kinematic factors and the arguments of the statistical functions that 
applies to the coUinear fringes, where either |x| ^ M or \y\ ^ M as well as close to 
the point of the soft divergence, where both \y\ ^ M. Using these expansions, in 
Section [531 we demonstrate the following points: 

• On the collinear fringes, where x = or y = 0, the total integrand takes finite 
values and is hence integrable, as it is expressed by Eqs. (IllSp and (11190 . 

• The point x = y = requires special care, since the soft gauge-boson singularity 
coincides with the Bose divergence of the gauge-boson distribution function. We 
show that the integrand behaves 1/ for {x, y) — )■ 0, such that this isolated 
singularity is integrable as well. A technical point is here that contributions from 
different regions of the angular integrations must be averaged. 

• The divergence from the Higgs boson distribution functions is integrable in the 
principal value sense. 
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In order to show the cancellation of the singularities on the coUinear fringes, we need to 
assume that Higgs bosons, leptons and gauge bosons are in thermal equilibrium. As a 
consequence, the coUinear sphtting processes are in equilibrium as well, what leads to the 
detailed balance relations f lll7p . that are essential in order to show that the integrand is 
finite on the fringes. 



3 NLO Neutrino Production Rate 



In this Section, we present the final results for the neutrino relaxation rate in terms of 
manifestly IR- and UV-finite integrals. The explicit expressions for the various terms 
are given in the subsequent Sections. The total relaxation rate in the plasma frame is 
[cf. Eq. 



2p' 



nO 



tr 



1 
2f 



tr 



1 LO 



+ tr 



1 wv 



tr 



1 VERT 



(19) 



Expressions for the leading order term 
function-type contributions are given bjo 



LO 



can be found in Ref. [2^]. The wave- 



tr 



wv 



^^wv,vac _j_ ^sca,vac^ _j_ ^wv,T^O _|_ ^sca,T7^0 (20) 

I ^^vac,col,wv _j_ ^vac.col.sca^ _|_ ^^HTL,col,wv _j_ ^HTL,col,sca^ 
_|_ jTvac.fin _|_ jrHTL.fin _|_ jrT^O 



and where the B terms correspond to radiation from the Higgs boson and J-" to radiation 
from the fermion. The various wave function, scattering, vacuum, coUinear and HTL 
(sub) contributions are isolated and collected in a way that the IR divergences cancel 
within each of the parentheses. The UV divergences cancel among the vacuum parts 
against the counter terms. The expressions for the various B and J-" contributions are 
given in Section HI 



^Due to Lorentz covariance: tr 



tr 
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The vertex-type contribution to the relaxation rate is given by 

27r 1 oo oo 

/^^(p-«)tr[^^^(p)]™ = -^^ I d{^-i^) I dcosg j dxj dyj^^\x,y) 

-10 

1 

-1 

1 

- —Y6Y / dcosg fe{E+)U{E_), (21) 

-1 

with 

J^^\x, y) = \ ( y) + r^\-x, -y) + r^\-x, y) + r^\x, -y) (22) 

+ JiT (^^ 1/) + y) + (^^ y) + ^2-''*(a^, 2/) 

where the J^^'^^{x,y) terms correspond to 2 -H- 2 scatterings and 1 ■v^ 3 decay and 
inverse decay processes with the kinematic region^ indicated in Figure [6l and J'^^^{x, y), 
J2±^{x, y) and J^^^{x, y) correspond to vertex corrections with an on-shell gauge boson, 
Higgs boson and lepton, respectively. The vacuum- type contributions J^^^'^^'^{x,y) and 
jj-vort.vacj^^^ ?/) are subtracted within J^^^{x, y) to make it UV finite (in addition to being 
manifestly IR finite), and they are cancelled by the vertex counter-term 5Y . The Fermi- 
Dirac distribution in the front of Eq. (1211) results from using the KMS relation (1A4I) . 
The explicit expressions for the integrands are presented in Section [51 In addition to the 
kinematic x, ^/-variables, the integrands depend non-trivially also on the angles — ip 
and p depicting the orientation between loop momenta and the plasma vector u, and E± 
are given by Eq. ( l82i) . 

4 Wave-Function Contributions 
4.1 Radiation from the Higgs Boson 

We now present the technical details that are omitted in Section 12.21 where a more 
qualitative overview of the present approach to the regulation and cancellation of the IR 
divergences is provided. 

^We have recast the integration area to positive x,y > 0, hence the exphcit minus signs in the 
arguments of contributions corresponding to regions II and III. 



18 



For a massless scalar radiating a gauge boson with the IR- regulating mass A, the 
vacuum and finite-temperature contributions to the hermitian self-energy are 



-rjH,va.c/- N G f d^q f 



(27" 



-(27,)^S\p~k-q){p+kY 
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(23b) 
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+4G 



(2^ 



27r5(fc2 _ a2)/^(A;0 



where G as defined by Eq. ([2]) encompasses the gauge coupling constants and u is the 



plasma vector defined in Eq. ([6]). The first of the integrals for n^'™'^ and 11 
respond to the sunset diagrams, the second of the integrals to the seagull diagrams (c/. 
Figure ED. 

A possible way of evaluating n^'™'^ is to go to the frame where p = and to introduce 
a momentum cutoff A. Then, 



cor- 



9(p' 
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G 
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G_ 
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log - + log2 + - 



A 



( log - + Ci ) + C2+p'6Z^ + 6mi 



(24a) 
(24b) 



where and Sm"^ are field-strength and mass counterterms. Alternatively, of course, 
an effectively equivalent result may be obtained using other regularisation procedures, 
e.g. dimensional regularisation in 4 — e dimensions, such that the leading terms in are 



Uf-{p) 



_G_ 

4^ 



fi 1 



1 



p' logY + T^Ae + - +p^6Z, 



(25) 



where 



A. 



7e + log(47r) 



(26) 



Notice that, as it is well known, the seagull graph vanishes when evaluated using dimen- 
sional regularisation. For both regularisation procedures, the UV divergences should 
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be cancelled by the field-strength renormalisation SZ^p. In the following, we effectively 
account for this by replacing A — )■ A, which takes a finite, renormalisation-scheme de- 
pendent value. 

For the part that vanishes as we take T — 0, we set A = and obtain 



n 



rp2 

(p) =G — a 



47r2|p| 



X 



rf|k| 



log 



2|k|/ + 2|k||p| 



p2-2|k|p0-2|k||p| 



log 



2|k|/ + 2|k||p| 



p2 + 2|k|p0 - 2|k||p| 



(27) 

/B(|k|). 



This is IR finite, such that it is indeed justified to take A = for this contribution. 

For the spectral self-energy, we perform the split into vacuum and finite-temperature 
parts accordingly. Notice that due to the CTP Feynman-rules, the spectral self-energy 
only receives sunset and no seagull contributions. The zero-temperature part is 
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{k) 



G- 



,2A;2- A2 k^ 



IGtt 



— ^t9(A;2 - A')signfc° ^ -^-^^^(A;^ - \^)signk' , (28) 
A;2 Stt 



where for the approximation, we neglect numerator terms of order A^. The Heaviside 
f^-function occurs due to the mass threshold of the would-be process of a scalar boson of 
mass square fc^ decaying into a massless scalar boson and a gauge boson of mass A. (The 
minus sign in front of the rate is not problematic, because an on-shell scalar particle 
cannot actually change its mass by radiating a gauge boson when the gauge symmetry 
is unbroken.) For the finite temperature contribution, we obtain 
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for P >0 



for A;2 < . 



(29a) 



(29b) 



We aim to calculate the correction B [defined in Eq. f lT^ ] by gauge-boson emission 
from the Higgs boson to the relaxation rate of right-handed neutrinos. The term that 
depends on 11^ can be identified as the wave-function correction (superscript wv), the 
term depending on 11^ as the correction from 2 -H- 2 scatterings and 1 -H- 3 decays and 
inverse decays, that we collectively denote as scatterings (superscript sea), such that we 
may write 



i3 = + B' 



(30) 



For the following calculations, we choose M = p° > 0, for definiteness, and note that B 
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should be odd in p^. The wave-function correction is given by 



(31) 



fcO = -pO±|k| 



X 



d 



-\ n?(fc) [f,[[p + k)-u) + U{k-u)] 



where we have made use of Eq. fll3al) . Clearly, we obtain only contributions for ± = 
±' = + when > 0. According to our decomposition above, this can be written as a 
sum of a term depending on n^'™'^ and 11^'-^^'', 



^wv,vac _|_ ^wv ,T^0 



(32) 



The contribution B^^'^^^ can be evaluated numerically, while g'"^'™'^ inherits the IR 
divergence from 11 J'™^, Eq. f l24bp . We can cast it to the form (|p| = 0) 



g^Y^GM'^ , A 



2,^4 log - / dn [feiip + k)-u) + Uik ■ u)] . 



(33) 



The scattering correction is 
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(27r)3 ^ (Pj^ 
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n^(*) [M(p + k) ■ t.) + Mk ■ (.)] . (34) 



Following the splitting of U-^ into vacuum and T 7^ contributions, we decompose as 
well 



gsca gsca.vac _j_ ^sca,TjtO 



(35) 



First, we show that g^^^a.T^o yigj^jg ^ finite result when A — )■ 0. For this purpose, note 
that for small — |k||, one may expand for k"^ > 
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(36a) 



and for k"^ < 



2k^--^lo, 



k|-A;0\ \k\-k^ 



+ 



(36b) 



Therefore, within ^'^ca.T^o^ ^g^^^ integrate over the singularity at /c^ = in the principal 
value sense. 
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Next, we demonstrate that jS^^^'"'^^ depends on A in a way that cancels the A depen- 
dence within ^wv'^c. For this purpose, we evaluate (p = 0) 



gsca,vac ^ _ g^y2 ^T^^^^^P"' ^ 2/|k| - A2)signfc° (37) 

\ J 1,0— „0-LI1,l r \ II 
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While for ^'^ca.vac ^ set A = and perform the integral numerically, we proceed 

with as 

|(m-a2/m) 

^sca,vac+ ^^^^2 ^ ^ ^ ^^^'j'j , j,^, [MiP + k) ■ u) + U{k ■ u)] (38) 



8tt J (27r) 

iog(M - 2|k|) (/,((p + A;) ■ «) + u{k . u))] 

+ J d\k\\ogiM-2\k\)^k^feiip + k)-u) + Uik-u)) 



Therefore, the terms oc logA^ cancel between and B^'^^. We emphasise that this 
cancellation works out, because n^'™'^ and n^'™'^ are imaginary and real part of the 
same analytic self-energy, which implies relation (fT6|) . 

While these expressions are already defined within the text above, we finally collect 
the explicit expressions for the various IR-finite contributions: 

/d^k 1 O 

j2^^MTm^ ^^'^^"^ ^'"^"'^^ ^'^^ ■ ''^^ ' ^''^ 

where k^ = —M — Ikl , 



s--^^°=,.y^|^ Y: ^,ii^'^^'{k)mp + k).u) + u{k.u)] (40) 



d^k ^ 2M 

A;0 = _A/±|k| 

with n-^''^^° given by Eqs. ([29]) and 

Bw..-^o =<,^,K= / (^2rf(M - 2|k|)^ (l + l-ij) [/,((p + i) . .) + Mk . .)] , 

(41) 

where fc" = —M + |k| = —M/2. In the last expression, we have substituted Il^''^^^{k), 
Eq. fl7r|) and have made use of the fact that the non-HTL part of n^'^^°(/c) [i.e. the 
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integral term in Eq. fj27|) ] and its derivative with respect to vanish for fc^ = 0. We 
also note that corresponds to the phase-space suppression due to the asymptotic 

thermal Higgs boson mass (note the minus sign from the statistical factors in the square 
brackets, as /c° < 0). Together with the IR-divergent pieces, these expressions above can 
be used in order to calculate the the NLO wave-function corrections to the right-handed 
neutrino production rate summarised in Eq. (!20|) . 



4.2 Radiation from the Fermion 



We now adapt the same strategy as for radiation from the scalar propagator to radiation 
from the fermionic one. As a complication, besides the cancellation between scattering 
and wave-function corrections from vacuum loops, there is also a cancellation between 
scattering and wave-function hard thermal loop (HTL) contributions. 

We express the leptonic self-energies as = PjCi^T^i^. In the approximation of 
massless particles in the loop, the spectral self-energy is given by 
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The bar on L[^^(y^,y) indicates that from this finite temperature term, the HTL con- 
tribution, which is purely thermal as well, is subtracted. 



23 



To this end, we also need the HTL-type contribution to the hermitian self-energy 



^i?,HTLj 



and the vacuum contribution 





lkl 




lkl 



16 k 



log 



a 

8^ 



A;0- 


h |k| 


A;0 - 


- |k| 




A 









GT^k' 



or, in dimensional regularisation. 
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The term SZ^ is a field-strength renormalisation that cancels the divergences as A — )■ oo 
or e — )■ 0, but is scheme dependent otherwise. In the following, we again implement its 
effect by replacing A — )■ A, where A is finite. The non-HTL contribution to the finite- 
temperature hermitian self-energy jg finite and can be expressed in terms of 



a one- dimensional integral, a classic result that can be found in Ref. [39 
The leading correction to the fermionic spectral function is 
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and the corrections to the Wightman functions are 
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In the following, we attach the same superscripts as for the functions l\^2 to various 
quantities, in order to indicate which loop terms they originate from. 
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The correction to the neutrino relaxation rate follows from Eq. f l9bp 
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where the particular contributions to J-" are defined below. 
The Dirac trace is evaluated as 
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As a consequence of Lorentz invariance, 

^AH.vacMj^^) oc fc/^, which may be verified 
explicitly by inspection of Eqs. fl42ll43ll45p . Therefore, we define 
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We have included here the threshold condition form the gauge-boson mass A, that is 
needed for the infrared regularisation. We obtain 



-4(p2 + e)k ■ S/'^'""^(A;) + A^p ■ S/'^'™=(A;) 



-2(A;V + (fc2)2)E/'^'™^(A;) . (53) 



The term oc {k'^Y gives a finite contribution to J^™c,fin ^^^^ evaluated straightfor- 
wardly. The term oc p^k'^ gives rise to jrvac.coi ^ jrvac,coi,wv^jrvac,coi,sca_ r^^^ wave-function 
contribution j^ac,coi,wv originates from S^'™^ and the scattering correction j:'vac,coi,sca 
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from S^'™'^. These terms individually contain coUinear divergences that are regulated 
by the gauge-boson mass A and that cancel when added to j^™c,coi^ demonstrate 
now. We find 
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where = p° — |k|. The scattering corrections yield 
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The logarithmic dependence on A is therefore cancelled cancelled with j^vac,coi,wv^ -g^ f[54|) . 
This is again a consequence of the fact that ^j^g i^j^e anti hermitian and hermitian 

parts of the same analytic self energy, evaluated at the two-particle branch cut. 

The HTL-type contributions are decomposed into terms originating from t^f and t^f 
as well, J^HTL _ jrHTL,wv^jrHTL,sca^ observe that A; ■ E'^'«TL(fc) = and ■ S^'^^tl(A;) 
is finite for k"^ = 0, whereas p ■ S^'^"'"^(fc) has a logarithmic divergence for fc^ — )■ 0. 
The terms that are oc p ■ S^'^"'"^(/c) and oc p ■ T,'^'^'^^{k) in Eq. fl50|) therefore lead 
to coUinear ly divergent contributions to J-"^"""^, and we denote these by j^htl,co1,wv 
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jrHTL,coi,sca^ whereas the terms oc k ■ T^^'^'^^^k) and oc k ■ S'^'^"'"^(/c) are finite and referred 

to as J^HTL,fin,wv jrHTL,fin,sca^ g^^j^ ^^lai 

jrHTL,col _jrHTL,col,wv _|_ jrHTL,col,sca (57a) 
jrHTL.fin _j7HTL,fin,wv _|_ jrHTL.fin.sca (57b) 

Tfie finite terms can be integrated in a straightforward manner, while the collinearly 
divergent contributions must be regularised, with a cancellation of the dependence on 
the regulator in the sum of wave-function and scattering contributions. In order to 
regulate these terms, we set the gauge-boson mass A = 0, but we introduce a lepton 
mass rrii. Then, we find that 

^HTL,coi,wv ^ _ g^y2 I ^^27r5(/ - |k| - y/k^ + mj) (58) 
^ ^ k'p- E™(fc) [1 - Mk . u) + U{{p - k) . u)] 



4A;o2 dko 



and 



2 P 

log dn[i- h{k ■ u) + U{{p - k) ■ u)] 



^ - ^'"^ y (2vr)3 2|k|(P-m|)2^^ ^^^^ 

x{l-h{k-u) + U{{p-k)-u)] 



2^2 



dn< [log (m^ + (2|k| -M)M) 



oo 

/ 

M 
2 



X [I - f,(k ■ u) + Uiip - k) ■ u)]]^l^-J^ 
d\k\ log (m,2 + (2|k| - M)M) ^ [1 - /^(^ ■ + /^((P - k) ■ u)] 



Notice that there is no + contribution from the sum over ± in Eq. because E^^^-^{k) 
vanishes for positive k^. The logarithmic dependence on mj therefore cancels when 
adding scattering and wave-function corrections, Eqs. (!59l) and (!58|) . Once more, the 
prefactors for the cancellation match because T,^"^^-^'^ (k) are anti hermitian and her- 
mitian part of an analytic self-energy evaluated at the branch cut. Notice, that the 
HTL-type contributions diverge logarithmically when — )■ 0. In that situation, one 
should use the resummed lepton propagator [l|-y, [351 . 



The result for J^"^^^ [that is defined as the contribution to J-' from and from 

fyA,T^o ^ which in turn results from replacing Jq — )■ 1^"^° and Ii — )■ 7^"^° in Eq. fH2]) ] can 
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be obtained by integrating over the pole at A;^ = in the principal value sense, what can 
be verified by checking the limiting behaviour of lQ^^{y^, y) and l[^'^{y^, y) for y^ — )■ y, 
i.e. that these functions are continuous at that point. 

For completeness, we again list the explicit expressions for the various IR-finite terms: 



vac, fin, wv 







(60) 



-J y ' \ \ fcO=M±|k| 

x[l-h{k-u) + U{{p-k)-u)] , 



(61) 



' (27r)^2|k| M + 2|k| Stt ^ ■'^^ i^Ji'Wi' ) )\ 



fcO=M+[k| ' 

(62) 



jrHTL,fin,v 



d^h 1 

(2;0^M2vr5(M-2|k|) 
GT^ 



Id 1 
+ 



dk^ 4k^^ 



X (^M^ + - k^) ^ [1 - Mk ■ u) + /^((p - k) . u)] 



where /c° = 



^HTLjfin.sca g 



(63) 



(64) 



d^k 1 
(27r)3 2|k| 



2^2 



2{p' + e)k ■ T.f^^\k) - 2ep ■ T.t''^\k) 



feO=^/_|_|k| 

x[l-f,{k-u) + U{{p-k)-u)] 
+27r5(M-2|k|) ' ^ 



(65) 



2k^^ dk^ 2A;03 



^) ((/ + ^^)A:-Sr^°(fc)-A:^p.Sf-^° 



(A:) 



x[l-/,(A;-M) + /^((p-fc)-M)] 



1.0— Ai 
J K — 2 



where jr^^-fi" = jp'acfin.wv^jrvacfin.sca^ exprcsslous can be substituted into Eq. ([20|), 

in order to obtain the production rate for right-handed neutrinos. Notice also that the 
sign of j^HTL,fin,wv jg consistent with the phase-space suppression due to the asymptotic 
thermal lepton mass. 
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5 Vertex Contribution 



5.1 Scatterings or Real Emissions 



The 2^2 and 1 o 3 contributions (involving the real emission of gauge bosons, here 
collectively referred to as scatterings) to the thermal production rate of follow from 
the first term of Eq. flTSl) , 



tr[^i^< (j9)]^^^ = Ag^GY^ 



[2(A;-j9) + g]''tr 



(27r)4 (27r)4 ' 

xiSf{kh''iS<{k + g)iA<(-A; + p)iA>,(g)iAj(-A; + p - q) . 
In order to bring the integrand to a symmetric form, we replace the momenta 



(66) 



q ^ —r 



p r 
^^2-^+2 



(67) 



which takes the effect 

, p r 

-k + p^- + q-- 



, P 

fc + fl— !■ q 

^2^2 



P 



-k + p-q^- + q + -. {6i 



Furthermore, the on- shell delta- functions of the Wightman propagators force 

P 



2+^-2 
p r 
2 ^ ^ ~ 2 



-2 =X\ 



where A again is a fictitious gauge-boson mass, such that this contribution reads 

d'^q d^r 



X ( p 

p 



(27r)4 (27r)4 



(69a) 

(69b) 
(69c) 

(70) 



+ 



+ q + - 
L2 ^ 2J 



p r 
- -q + - 
2 ^ 2 



x2n6 
xsign 



p r \ 2 

q 

^ 2 



2 

M 



p r 
2 ~ ^ ^ 2 



^0 



2n6 



sign 



PV 



p r \ 2 

- + q- - 
2^2 



M 
P 



+ q 



(f-^+i) (f+^+i)- 

iign (r°) 
n) [1 + fA{-r ■ u)] 
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We integrate first over dr^, dq^ and dcosi} using the three on-shell delta-functions and 
then change the integration variables |r| and |q| back to r° and i.e. the integration 
measure changes as 



d^q d^r 
(27r)4 (27r) 



-27r(5 



p r 
2-^-2 



(71) 



1 



2tv 



2-E 



(27r)5 

q0=±'|50| 



d(f d-ip dcosg / q^(i|q| / r^(i|r 



1 



-1 



2Vr2 + A2 4|g0| |q||r| 



27r 



27r 



1 1 



(27r)^ 



dip J dip J d cos ^ y (ig° y dr^ 

-1 X I 





^ sin "(9 sin ^ 






. 


^ sin g sin ^ 


r = r 


sin ■}) cos (y9 


> q= |q| 





P 


sin g cos t/' 




y cos -i? y 






y cos g J 



where g denotes the angle between q and u. More specifically, we choose to parametrise 
the angular dependences as 



(72) 



The integration domains X are best determined from distinguishing between the 
various kinematic situations, which are sketched in Figure [61 Besides, we take within 
these regions a positive integration measure, such that we do not obtain explicit minus 
signs from the Jacobians. The regions follow from the three on-shell constraints. Notice 
that these conditions may be combined to obtain 

, , 1 



|q| 



and when using this clS db relation for cos 

2 



r 



M 



-.0 



M 



_ _l_ gU J _|_ |q| (,Qg^ _|_ _ q- 



0. 



(73) 



(74) 



The relation — 1 < cost? < 1 is then one of the constraints that must be imposed in order 
to determine the regions of integration in Figure [61 Indeed, for A = this inequality 
reduces to: 



02 



4 



> 0. 



(75) 



Besides, we note that the on-shell constraints imply that the invariant momentum 
squares of the off-shell propagators are 



P 2 



T) r 
- -q+ - 
2 ^ 2 



--M (2g° + r°) , 
--M (-2g° + r°) 



(76a) 
(76b) 
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This suggests a change of variables to 



a; = g° + , (77a) 

y = -q' + \r\ (77b) 

which has the Jacobian —1. The integration regions in terms of these variables are 
indicated in Figure |6] as well. The collinear divergences are located on the fringes of the 
integration region where either x or y vanish, soft divergences are at the coincident point 
X = y = 0. 

We can therefore recast the scattering contributions to the neutrino production as 

2tt 1 

tr[^i^<(p)]-^=4^7^Gy2l_i_ J di^-i;) j dcosg j dx j dy (78) 

-1 Xy 

fM \ . [M \ ^^-2M^ + 4:Mx + 2My + 4xy 
xsign(x + yjsign — x \ sign — y \ PV- 



2 y V 2 y xy 

xfe {{p/2 - q - r/2) ■ u) U ((p/2 + q - r/2) ■u)[l + fA{-r ■ u)] . 

The statistical functions can be evaluated using above relations. In particular, the scalar 
products that appear in the arguments are given by 

p ■ u =p^ , (79a) 
(7 ■ M =-(x — — Iql-^ cos p, (79b) 

IpI 

r ■ u ={x + y)j^ — |a; + y\j^ [cos cost? + sin sin cos(v9 — ip)] . (79c) 
with |q| as in Eq. fl73|) and cost? given by Eq. flTijl . Close to the collinear edges, where 
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a; = OVy = 0, we can evaluate 



- - g - - j ■ " ~ ( 2 ~ jvf j + IpI ^ sign(M - x)sign(M - y)) (80a) 



+-^|p|sin^cos(v?-V), 
' 9 - o • " ~ o - 77 (P° - IpI cos sign(M - a;)sign(M - y)) (80b) 



2^2/ \2 M 



■Jxy 

+ -^|p|sin^cos(¥?-V), 
r ■ u ^ "I '^^^ ^ sign(M — x)sign(M — y) (80c) 



IpI sm gcos[(p — ip) 



M 

Note that Eqs. f lHOj) are understood as leading order expansions and hence the ~ ^Jxy 
corrections are applicable only when both x ~ y ~ 0. As we aim for rearranging the 
various contributions to the vertex-type correction into a manifestly finite integral, we 
have dropped here the dependence on the regulating gauge-boson mass A. We have also 
let cos Q — )■ sign(M — x)sign(M — y) cos q for later convenience. 

In order to arrange for an IR-finite integral, where the collinear divergences on the 
fringes a; = and y = Q cancel, it is useful to define the integrand 

J'-^x, y) =K''\x, y)h ((p/2 - q - r/2) ■ u) U ((p/2 + q - r/2) ■ u) (81) 
X [1 + fAi-r ■u)]S{x,y) 



|a;|<MV|y|<gM \\ M J J \\ M 

where s = sign(M — x)sign(M — y), 

E± = ^{po±\p\cosg), (82) 

and S{x,y) is a step function defining the support of the integral in accordance with 
Figure |6l 

S{x, y) =^{x)^{y)^{M/2 - x - y) + ^{-x)d{-y) (83) 
+ {^{x)^{-y) + ^{-x)^{y))^ {x + y- M/2) . 
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For later convenience, we have separated the kinematic part in the factor 



r^sca/ N . / \ ■ (M \^^-2M^ + 4Mx + 2My + 4xy 
K [x, y) =sign(a; + yjsign 1 y - a; 1 sign 1 y - Z/ 1 PV . 

(84) 

In the approximate form for J^^^, we have made the restriction to |x| <^ M V \y\ <^ M, 
where V denotes the exclusive or, which justifies dropping the ~ ^Jxy corrections. We are 
going to use this expression in order to show the cancellation of collinear divergences. For 
the soft divergence at x = ?/ = 0, there occurs a complication due to the Bose divergence 
of the gauge boson distribution, and we will provide a separate discussion. 



5.2 Vertex Correction with on-Shell Gauge Boson 

We now bring the vertex contribution with an on-shell gauge boson to a form that can be 
matched with the contributions from scatterings. It is useful to reparametrise in Eq. f[T5|) 



such that the second term in Eq. f|T8|) turns into 



tr[fit^ip)]r'={3gl + gl)Y 



[-k - qYti 



(85) 



(86) 



(27r)4 (27r)4 

xiSfip - q)Y^S<{p - k)iA<{k)iA^ {qyKM - k) 

We perform the integrations over dk^ and dq^ making use of the on-shell delta-functions. 
When we define 



fc± = (±|k|,k), gi± = (A:°±u;,q) 



v/(q-k)2 + A2, 



we obtain 
tr[|ii^ 



2-K 



< /^M vert 
1 



'^^"^^'16(^5^ / j d cos gjdu j d|q|M 



(87) 



(88) 



xPV 





2M^ 



^ Qmin 



[(M/2 T w)2 - q2] [(M/2 ±00)- q^] 
2M2 M2 



(M/2 ^ w)2 - q2 (M/2 ± w)2 - q2 _ 
X [1 + 2/A(|(gi± - k+) ■ u\)] mp - k+) ■ u)U{k+ ■ u) , 



where 



imax,mm 



2 



(89) 
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Here, we parametrise the angular dependences as 



q= |q| 



( sin -i? sin (/9 ^ 

sin 1? cos if 
\ cos -d J 



k= Ikl 



/0\ 





u = 



M 



( 



sm Q sm 



sin Q cos -0 

y cos Q j 



(90) 



In order to evaluate the distribution functions, we use 
k+ ■ u = -{po — |p| cos g) — , 



{p-k+)-u^ -{po + IpI cos g)^E+, 



qi±-u^Po ( - |p|^(cos^?cos^ + sin^?sin^cos((^- V')) , 



COS'J? 



M\q\ 

Next, we change variables to 



(91a) 
(91b) 
(91c) 

(91d) 



X 



2M 



M 



-uj - q 



y 



2M 



M 



+ 00 -q^ 



(92) 



Notice that ,x = corresponds to the pole of the lepton (Higgs) propagator for ± — > + 
(d= — > — ) and that vice versa, y = corresponds to the pole of the lepton (Higgs) 
propagator for ± — )■ — (± — )■ +). The inverse transformation gives 



|q| = ^ + AM{x -y) + A{x + yf , 
u ^ x + y , 



(93a) 
(93b) 



and the Jacobian is — M/|q|. We express the vertex contributions as (we take A — >■ at 
this point): 



27r 



-1 /» n /"OO noo 

tT\^it<{p)]J-'^^^4g^GY'^^^Jd{^-i^)Jdcosgj^ dy dx (94) 



X 



where for the present terms 



J7r(^,2/) =KTi^:y) [1 + 2/A(|(gi± - k+) ■ u\)] fe{E+)U{E_) , 



(95) 
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with 



+ 2Mx - My ^, -2My + Mx 



xy 



xy 



(96) 



In this form, we can directly compare with the terms from real emissions. When 
either or both, x, y are close to zero, we can approximate 



|q| ^ -|M + 2x- 2y\ 



(97) 



^ _j_ y ^ y ^f%^ 

(?i± -h^)-u'^ iPo^TT IpI cos^— 2— — IpI sin^cos((y9 - ^) . (98) 



M 



M 



M 



Expanding the arguments of the distribution functions close to the surfaces x 
Vy = 0, we obtain 



Jir{x,y)^Kri\x,y) 



2x 



2y 



1 + 2/a —E^ + -i-K 



M 



M 



ME^)U{E_) . (99) 



5.3 Vertex Correction with on-Shell Higgs Boson 

We now consider the third term of Eq. flTSl) . which using the reparametrisation 
becomes 



J (27r)4(27r)4 
xiSfip - q)Y^S<{p - A:)iA<(A:)iAj(g)iA2.(g - k) 

Making use of the on-shell delta-functions, we find 



(100) 



tT[0^M]r' =^^-^^'^^(L)^ ^ C '^^^ -^)Jd^osQjdu I dq^ (101) 



xPV 



(M ± a2q){M/2 ± aq + uj){M/2 ± aq - oo) 
AM 2M^ 



M±a2q {M/2 ± aq + uj){M/2 ± aq - u) _ 
X [1 + 2/^(1 (g2±a - k+) ■ u\)] feiip - k+) ■ u)f^{k^ ■ u) 



where we define 



M 

a = sign I Y - 1/ 



(102) 
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Note that in the expression f llOip . g is a real number, not a four- vector. In terms of the 
original integration variables, the new variables are given here by 



v/^^+P, g= v/(k-q)2. (103) 



We again change variables according to Eqs. f l92|) and with |q| — and decompose 
the present contribution as in Eqs. f lM|) . Notice that the IR divergences only occur for 
the contributions where ± — )■ — , such that apart from the Bose divergence in the Higgs 
boson distribution function tr[|ii|^^(p)]™''* is finite. The integrand for Eq. flM|) with i = 2 
then is 

Jr±\x,y) =Klf{x,y) [1 + 2/^(|(g2±. - A;+) ■ u\)] , (104) 

where /c+ ■ u and {p — kj^) ■ u are again given by Eq. ( l9T|l and 



Klf{x,y) =-PV 



+T7^^ (105) 



(M ± a2q) (M/2 ± aq + uj) (M/2 ±aq-Lo) M± (j2q 
2M^ 



(M/2 ± org + uj){M/2 ± aq - u) 

We do not explicitly substitute x and ?/ here for q and w as this would give rise to a 
lengthy expression in this case. 

In order to compare with the other contributions, we notice that when either or both, 
X and y are small, we may approximate 

— q + u] ( q — uj) ~ — 4x?/ , (106a) 



2 / V 2 

M - 2g ^2{y - x) . (106b) 

The relations for constructing the arguments of the distribution functions are 



(1 q \ y/cu'^ — 
- ± — j - IpI — (cos cos ^ + sin sin g cos{ip -'?/')), (107a) 

cosi? = \ — . (107b) 

When either x, ?/ or both are close to zero, we may approximate (again, we take A — ?■ 
here) 



{q2± - k+) - u^ {^ + ^^jf-j (±Po + IpI cos ^) - 2-^|p|sin^cos((/? -?/^) . (108) 



36 



Using these approximations, we find for x <^ M y y <^ M 

AM{x + y) 



Tv'vertjappr 
^2+ ' 



Tv'vertjappr 
^2- 



x,y) =PV 
x,y) =PV 



(M + 2x - 2yY ' 

M{M + x-y) 
xy 



AM 



M + 2x-2y 



sign(x — y) 



J,Ti^, y) ^Klf^'-^'^ix, y)sign(M - 2y) 



(109a) 
(109b) 
(109c) 



X 



1 + 2/0 



1 + 



2x-2y 



M 



5.4 Vertex Correction with on-Shell Lepton 

Tfie tliird term of Eq. ( IT8l) in tfie reparametrisation ( l85l) is 

,2 f d'^k d^q 



tr hKii) 



'N 



vert 
3 



k - qYti 



(27r)4 (27r)4 

xi^f (p - qh'iS^ip - k)iA<{k)iAf (g)iA2.(g - k) 



(110) 



The kinematic situation is essentially the same as for the on-shell Higgs boson, which is 
why the relations fll06lll07lll08p apply to the present contribution as well. The integrand 
can be expressed as 



J^l'\x,y) =Kr\x,y) [1 - 2/,(|(g2±. - fc+) ■ n|)] ME_)U{E^) , 



where 



K;^^{x,y) =-FV 



AM^ 



+ 



AM 



(M ± a2q) {M/2 ± aq + uj) (M/2 ±aq-Lo) M± a2q 
AM^ 



[111) 



[112) 



{M/2 ±aq + u){M/2 ± aq - u) 



The arguments for the distribution functions can be constructed using Eqs. (1107^ . with 
the important exception that we replace cos — )■ — cos g in order to interchange <f->- 
E-. This change of variable leaves the value of the integral unaffected, but it is helpful 
in order to to combine the integrands to a manifestly IR-convergent expression for the 
right-handed neutrino production rate. 

The portion of the integrand that leads to IR divergences is approximated for x <C 
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M V 2/ < M as 



Tv'vertjappr 



^^vert,appr | 



{x,y) =0, 
x,y) =PV 



M(M + 2x- 2y) 
xy 



AM 



M + 2x-2y 



sign(a; — y) 



J,Tix. y) ^Kl^^'^'^'^x, y)sign(M - 2y) 



(113) 
(114) 

(115) 



1-2/, 



1 + 



2x - 2y 



M 



5.5 Infrared Finite Integral 



The remaining task is to combine the various pieces to an integral, that can be evalu- 
ated with a finite result (in the principal value sense) for each subset of the complete 
integration region. This is achieved by the combination 

2ir 1 



oo oo 



tr[^i^;(p)]^^^^ =^9u,GY' 



32(27r) 



d{v-^lj) / rfcos^ / dx / dyj^°^''\x,y) , (116a 



-1 







J'°'^\x, y) =i- 1 r^'ix, y) + r'\-x, y) + r''\x, -y) + jrJ\x, y) + jr\x, y) 



+:r^^^(x, y) + J^'^^x, -y) + J^'^^-x, y) + J^f{y, x) + j;i^\y, x) 
+m-y, -x) + Jr+\x, y) + J^l^\x, y) 
+r^^{-y, -x) + jrrHl/, x) + J^l^\y, x) 



+ + 2J^f (x, y) + 2J^f{y, x) + x ^ y \ + ^ ^ ^ + -k . 



(116b) 



(Recall that the support of J'^'^^ is given by Eq. fl83|) and shown in Figure [61) This 
expression may be integrated without introducing an infrared regulator. 

In order to show that this is integrable on the fringes where x — )■ V|/ — ?■ 0, we make 
use of the detailed balance relations applicable to coUinear splittings 

U{E) + U{E)U{aE) + U{E)fA{{l - a)E) =U{aE)U{l - a)E) , (117a) 

fe{E) - fe{E)f,{aE) + ME)U{1 - a)E) =fe{aE)fA{{l - a)E) . (117b) 

It then follows that for x — )■ (divergences in the Higgs boson [lepton] propagator for 
Eq. ffTTSil) [Eq. f lTT8b|) ]) 

xy {J^'^ix, y) + J'-%-x, y) + J^l^'^x, y) + J^l'\x, y)) ^ , (118a) 
xy ( J^^"(-l/, -x) + ^) + J^-\y. ^)) ^ (118b) 
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and that for y — > (divergences in the lepton [Higgs boson] propagator for Eq. flllQal) 
[Eq. ffTT9b|) ]) 

xy {r'\x, y) + r'\x, -y) + T,f{y, x) + j;i'\y, x)) ^ , (119a) 

xy{r'\-y,-x) + Jlf{x,y) + jr\x,y)) ^0. (119b) 

The sum J^^^{x,y) + J^'^^{—x,y) in Eq. flllSap makes it sure that the cancellation with 
vertex contributions takes place on the whole coUinear edge x = (for both y < M/2 
and y > M/2), and similarly in Eq. f lll9ap for y = in accordance with the kinematic 
regions, cf. Figured In order to obtain the manifestly IR-finite result for the vertex- type 
part of the neutrino production rate flll6al) . we need to add these contributions to those 
terms, that are integrable from the outset and finally symmetrise x ^ y, in order to 
make sure that the sum of all divergent terms cancels on both edges (x = and y = 0). 
Furthermore, we have symmetrised over the polar angles. 

The criteria (111811119 1 ) are in general not sufficient in order to show the integrability 
at the point x = y = 0, because the detailed balance relations only hold on the collinear 
fringes where x = 0\/ y = 0. We must therefore show in addition that x?/j7'*°*'^^(x, y) van- 
ishes on each trajectory for [x, y) — )■ 0. The problematic terms are the Bose divergences 
in J'^^'^ and JT^"^*- When using Eqs. fl81|95l) together with the approximations for the 
arguments of the distribution functions (I80II98I) . we notice the cancellation of the diver- 
gences in xy{J'^^'^{x, y) + Ji^^{x, y)). There is however a finite remainder that moreover 
depends on the direction of approach of {x,y) toward (0,0). This finite contribution 
follows from the terms that are linear in x, y and ^Jxy in the arguments of the distribu- 
tion functions (I80ll98p . For this finite term to cancel as well, we must add X^^^{—x, —y), 
symmetrise over x and y (as it is already necessary for the cancellation of the collinear 
divergences) and moreover, add the contributions from positive and negative values of 
cos((y9 — In summary, it follows that 

lim xyj''"''\x,y) = 0, (120) 

such that J'^°^^^{x,y) grows more slowly than l/(x^ + ?/^); and therefore it is integrable 
at X = ?/ = 0. In fact, we have verified by numerical evaluation that J7'*°*^'(x, y) ~ 
1 / \J x^ ^ y'^ for (x, y) — )■ 0. 

Similarly, care must be taken about the points (x, ?/) = (0,M/2) and (x, y) = 
(M/2,0), where the Bose divergence of the Higgs scalar coincides with the collinear 
fringes. Suppose we are close to the fringe where x — )■ 0. Inspecting the approxima- 
tion fll09p for J2±^ 1 we notice that the y integration can be performed oX y = M/2 in 
the principal value sense. The same is true for the collinear ly divergent terms oc 1/x 
in J7''^'^'^( sign (M/2 — y)x,y), Eq. ( IHTll . Therefore, we can perform the dAy integral over 
jrtotai(^^^ M/2 - Ay) + J*°*^'(x, M/2 + Ay) as well {Ay > 0). Due to the detailed balance 
relations flll8|119p . x times this integral vanishes for x — )■ 0, which is why the remaining 
X integration yields a finite result as well. We have also checked numerically that for 
(x,i/)^(0,M/2) 

x{y-M/2)j'°'^\x,y)^ 0, (121) 
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and moreover that J7'*°*^'(x, y) ~ + — M/2Y irrespective of the hmiting di- 

rection, which further confirms the integrabihty at (0,M/2) and (M/2,0) hy x ^ y 
symmetry. Furthermore, we have verified that the Bose divergences in J2±^ 
hne y = M/2 + x away from fringe x = are actually cancelled in that [y — M/2 — 
x)J'^"^^\x, y) — )■ and J'^°^^\x, y) — finite for y — )■ M/2 + x , and similarly for x y. 

5.6 Ultraviolet Regularisation 

The vertex-type correction flll6al) to the neutrino production contains UV divergences 
that need to be isolated and subtracted in order to obtain a remainder that is suitable for 
numerical integration. The UV divergences originate from the vertex contributions. In 
particular, the integrands fl95|104|llll) are not exponentially suppressed for large values 

01 X oi y because of the constant term in the square brackets. The contributions arising 
from this constant term are easily interpreted as the vacuum vertex-corrections to the 

2 — >■ 1 processes — )■ or £0^ — )■ A^, weighted by the statistical distributions for £ 
and 0. These UV-divergent contributions are analytically calculable, and are given by 
the the following sum of terms: 

J-t,vac(^^ ^) = ^ {Klf{x, y) + Klf{x, y) + iTJf (x, y)) h{E^)U{E.) . (122a) 

Now, we need to subtract this expression from the integrand that we aim to evaluate 
numerically, while adding terms that balance the collinear divergences that this subtrac- 
tion may induce. The vacuum-vertex corrections should only have IR divergences that 
cancel with those for real emissions in 1 — i- 3 processes. (Scatterings of the 2 -H- 2 type 
do not occur in the vacuum, in absence of incoming gauge bosons, leptons and Higgs 
bosons.) This implies that outside of Region I, J^^{x,y) should not yield coUinearly 
divergent contributions, which can be verified as 

x?/( J""''*'™"(x, y) + J^™'^*'™^(y, x)) = for (x = A y > M/2) V (y = A x > M/2) . 

(123) 

This is no longer true for Region I. In order to cancel collinear divergences there, we 
define 

jr™-(x, y) = K^'\x, y)^{x + y- M/2)f,{E+)U{E^) , (124) 

which we subtract in addition. This is just 1 — ?■ 3 rate weighted with the same statistical 
factors for leptons and Higgs bosons as J7'^'^''*'™'^(x, y), and it has the desired property 
that 

xy(J"^(x,?/) + j7™'"*'""'(x,?/) + xO?/) = forx = OV?/ = 0. (125) 
Now the integrand 

J^^\x, y) = J*°*^'(x, y)-^ [J^^^'^^^ix, y) + J™'^*'^^^(x, y) + x ^ y) (126) 
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is IR and UV finite and suitable for numerical evaluation. 

Finally, the subtracted terms must be calculated analytically and then added again. 
For this purpose, we note that the virtual T = corrections to the decay rate, evaluated 
with a cutoff u < A [from the parametrisations in Eqs. f l55|10ip ] are 



M2 



-21og^ 



A 



+ 



TT 



A 



A 



M 6 



+ log2 



(127) 



+YSY' 



Here, 6Y is a counter term that must remove the logarithmic divergence in A and the 
further details of which are determined by a particular renormalisation scheme. 

Alternatively, we can use dimensional regularisation in order to calculate the reduced 
vertex function 



V{p) =2g^GY'fi' 



Itr 



(27r)4- [{p - g)2 + ie][{q - kf - A^ + ie][q^ + ie] 



+ Y5Y2M^ 



128) 



--2g^GY'- 



8^ 



-2 log' + 3 log - + 2 log — 
M 6 fi fi 



A. - 



Y6Y2M^ 



In this form, SY must cancel the 1/e divergence and may otherwise again be renormalisation- 
scheme dependent. Integrating over the T = phase space, we obtain the zero-temperature 
spectral-function for the right-handed neutrino 



trk>Ei{p)r-J, 



"^'^ "^'"^ -i2n)^S\p -k- q)27rSik^)27rSiq')Vip) = ^ (129) 



(27r)4 (27r) 



2^713 

.M2 



-21og' — + — -31og- + 21og — 
M 6 fi ji 



A. 



57r 

1' 
2 



-Y5Y'- 



Note the factors from the two different diagrammatic contributions in the CTP approach 
(cut on the left or on the right). Comparing with the prefactors in Eq. (!94|l and using 
the definitions (lAll) . we can conclude that 



J' 



dx J dyj 





vert, vac , 



J 



vort,vac / 



x,y) = -2M^ 



A TT 



M 



-2 log^ — + — - 3 log - + 2 log A, - - 



M 6 



(130a) 
(130b) 



41 



The IR divergences are again isolated in terms involving log A, which cancel with corre- 
sponding terms in the vacuum vertex contribution fll30ap . 

Next, we need to obtain a matching expression for J^'^^'^^{x,y), which we have 
introduced in order to cancel the IR-divergent A-dependence. For this purpose, we 
evaluate the integral 



M M 
2 2 



dx I dy 



2L 2l£ 

2M ix 



M M ' 
-2 + 4— + 2— + 4 

xy y X 



2M^ 



^, 2 A vr^ ^, M 11 . 

-2 log^ — + — + 3 log 

^ M 6 ^ A 2 



:i3i) 



where the boundaries of integration derive from the condition —1 < cos'd < 1 with 
Eq. (O]). Hence, 



oo oo 



J^^'^'^'^^ = J dx j dyj^^'^'^^'^^x^y) 



(132) 




r2 



A 



M 11 



The final expression for the reduced self-energy is 



2tt 



oo oo 



tr[^i^< (p)]™ =4<7^Gy2_i_ J dicp-^) J dcosg J dx J dyj^^\x,y) (133) 

-10 

1 
-1 

1 

+— W J dcosg ME+)U{E^). 
-1 

Each of the explicit integrals in this expression is convergent, while the remaining UV 
and IR divergences are isolated within simple factors that cancel in the sum, as discussed 
above in this present Section. 

6 Discussion and Conclusions 

The results of this work show that the relaxation rate of right-handed neutrinos 
toward thermal equilibrium is finite to NLO in perturbation theory, i. e. when including 
the leading corrections from Standard Model gauge radiation. As we perform a fully 
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relativistic calculation, we generalise earlier results that are based on non-relativistic 
approximations js], ^ that are valid when M ^ T. In particular, the relative motion 
of the neutrino with respect to the plasma is kept general here and the full quantum 
statistical Bose-Einstein and Fermi-Dirac effects are accounted for. 

We have developed two somewhat different methods in order to handle the wave- 
function- and the vertex-type corrections. The treatment of the wave-function corrections 
is greatly facilitated, because well known and relatively simple analytical expressions for 
the thermal self-energies are available. We can isolate the IR divergences in terms of 
logarithmic dependences on the regulating gauge-boson mass A for both, the scattering 
contributions that rely on the spectral self-energy and the wave-function contributions 
originating from the hermitian self-energy. Eventually, we find that the IR divergences 
cancel [Eq. ([33]) with Eq. ([38]), Eq ([5l with Eq. §E> and Eq §E> with Eq. ([59])] because 
the hermitian and the spectral self-energies are the real and imaginary parts of the same 
analytic function evaluated at the two-particle branch cut. 

For the two-particle irreducible vertex-type corrections, the situation appears more 
complicated. Rather than isolating the IR divergences in analytic expressions for vari- 
ous contributions, we therefore manipulate these in such a manner that we obtain the 
integrand (]116b|) that is manifestly free from non-integrable IR and Bose divergences. 

Since the methods that we use in order to demonstrate the cancellation are very 
explicit and as we also discuss the subtraction of UV divergences, they readily suggest 
a method for analytic or numerical evaluation of the relaxation rates of right-handed 
neutrinos A^. In a future work, it would be of interest to verify the results of Refs. js], 0] 
by taking the non-relativistic limit of the present results and moreover, to perform a 
numerical evaluation of the production rate of without non-relativistic approximations. 

While the present work provides a method for evaluating the relaxation rate of right- 
handed neutrinos N, that derives from trhz*^^], cf. Eq. ([3]), the rates of non-equilibrium 



CP- violation from decays and inverse decays of N typically rely on different components 



of the self energy [2J, |29|, |31|-|33| . Another future task would therefore be the evaluation 
of CP-violating rates, in generalisation of the methods presented here. 

While we have focused in the present work on the production of right-handed neutri- 
nos, the basic topology of the self-energy diagrams that describe the decay and scattering 
rates and their NLO corrections is obviously shared with diagrams for other high-energy 
reactions in finite-temperature backgrounds. For that reason, it may be useful to for- 
mulate our method in terms of master integrals in such a way that the results could be 
more directly applied to other situations with different kinds of particles in the loops 
(different Lorentz tensor structures). One may investigate for example the processes that 
underlay the transport coefficients used in calculations for Electroweak Baryogenesis or 
corrections to the production rate of Dark Matter particles. It will therefore be interest- 
ing to explore the possibilities of applying the present methods to additional topics in 
Early Universe Cosmology. 
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Appendix 



We follow the conventions of Ref . 16| . The relations between the OTP two-point func- 
tions, (G stands for a propagator A or 5 or for a self-energy 11 or ^) are given by: 



= 


G'^ -G> = G< - G^ 


(advanced) , 


(Ala) 




aT_Q< = Q>_ qt 


(retarded) , 


(Alb) 




= i(G^ + G^) 


(Hermitian) , 


(Ale) 


G^ = ^(6-^ 


-G^) = -G<) 


(anti-Hermitian, spectral) , 


(Aid) 


G^ = -(G> 
2^ 


+ G<) = i(G^ + G^) 


(statistic) . 


(Ale) 



In terms of the functions bearing OTP indices, these can be expressed as 

G^ =G~^~^ , G^ = G (time ordered, anti time-ordered) , (A2a) 
G< =G+- , G> = G-+ (Wightman) . (A2b) 

The Wigner transform is defined by 

G(p, x) = j d'^re'P''G{x + r/2, x - r/2) . (A3) 

Under spatially homogeneous conditions, there is no dependence on the average spatial 
coordinate x and furthermore, we suppress the explicit average time coordinate t = x^. 
Equilibrium Green functions or self energies observe the Kubo-Martin-Schwinger (KMS) 
relation 

G>(p) = ±e^P"G<(p), (A4) 

where the plus sign applies to bosonic, the minus sign to fermionic two-point functions, 
and u is the plasma four-velocity. 
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The tree-level equilibrium propagators for the lepton doublet that we use here are 

iS<ip) = -27r5(p2)sign(/)PL/^/,(p ■ u) , (A5a) 
iS>{p) = 27r5(/)sign(/)PLj^^(l - • u)) , (A5b) 

iSjip) = Pl-^ - 27r6{p')PMi\P ■ «l) , (A5c) 



= -^^^^ - ^^Kp')PLm\p'\) , (A5d) 

where fi is the Fermi-Dirac distribution. The generalisation to non-equilibrium distribu- 
tions can be found in Ref. [2J]. We suppress explicit SU(2) indices, which are contracted 
with that of the Higgs doublet in the usual manner, using the two-dimensional anti- 
symmetric tensor. 

Similarly, we use the scalar equilibrium propagators 

iA<{p) = 27r5(p2)sign(/)^(p ■ u) , (A6a) 
iA>{p) = 27r5(p2)sign(/)(l + Uip ■ u)) , (A6b) 

iAj(p) = + 2n6{p')U{\p -ul), (A6c) 

^ p"^ + le 

iAj(p) = --^ + 2nS{p')U{\p -ul), (A6d) 

where is a Bose-Einstein distribution. The gauge field propagators are 

iA<,(p) = 2nS{p')sign{p'){-g,,)fAip ■ u) , (A7a) 
iA>M = 2n5{p')sign{p'){-g,,)il + /^(p ■ u)) , (A7b) 

iAj.(p) = + 27r6{p'){-g,^)U\p ■ n|) , (A7c) 

p^ + le 

i^Up) = + ^^Sip')i-9,.)fAi\p ■ n|) , (A7d) 

p"^ — le 

and Ja is a Bose-Einstein distribution. 

As a consequence of the Majorana condition, the self energy of the right-handed 
neutrinos inherits the property 

tNix,y)=CtUy,^)C\ (A8) 



where the transposition t acts here on both the CTP and the Dirac indices, cf. Ref. [2^ 
for a more detailed discussion that also includes the mixing of several right-handed 
neutrinos and the effect of non-zero chemical potentials for the leptons i. 
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